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Abstract. 

The n*'' derivative of a tensor valued function of a tensor is defined by a finite number 
of coefficients each with closed form expression. 



1. Introduction. We consider tensor functions on symmetric second order tensors, 
j3 ' S'ym— > Sym, defined by a scalar function f{x) of a single variable according to 



/(A) = ^/(a.)A„ (1.1) 

cn ■ 

where d is the eigen-index of A. The tensor A g Sym is arbitrary with spectral decom- 
position 

d f A ■ — ■ 

A = ^«,A„ A,A, = i ' J2^=^ (1-2) 

■ The specific case of tensors acting on 3-dimensional vectors, d < 3, is discussed in this 

! paper, although the results can be readily generalized. 

' Derivatives of /(A) are defined by the expansion 

/(A + X) =/(A) + V/(A)X + iv(2)/(A) : XX + ^V(3)/(A) : XXX 

+ . . . + -^V<"V(A) : XX . ..X + ... (1.3) 

n terms 

The n*'* derivative V("V(A) is a tensor of order 2{n + I) which contracts n-times with 
the second order tensor X to produce a second order tensor. The first derivative, or 
gradient, is 

V/(A) = J2 /i/^A,HA„ = <j (1.4) 
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where denotes the outer tensor product, defined in Section [2l The identity (|1.4[1 is well 
known and has appeared in various formats. The first explicit presentation I am aware 
of is due to Ogden [10] who defines a fourth order tensor C} — df{A)/dA (with slight 
change in notation). Ogden gives the components >C,]jj,; in terms of the eigenvectors 
of A, i.e. — C}jf.iai ® aj ® (g) a;. These coefficients are related to those of (|1.4p 

by the fact that = a^ (g) a^ (no sum) when d = 3 and hence Cj^^j^ = f^l\ Cjj^j — f^^\ 
The fundamental result (II. 4p was also derived by Carlson and Hoger [T], although the 
present notation is based on [T3] . 

Ogden fTIJl (Section 3.4) also presented the second derivative. In the present notation 
it is 

-V(2)/(A)= ^ /g^A,HA,KAfc, (1.5) 

i,j,fe=l 
fiii ^2"^ i^i)^ 

A2) ^ fi^j) - /(«») - («j - "»)f («0 

j(2) _ Ifjaj) - fiai)]{at + «j - ^ak) ~~ [/(«») + /(«j) - 2f{ak)]{aj - a») 

2{ai - aj){aj - ak){ak - ai) 

Ogden's expressions (|1.4|) and (jl.Sp are special cases of the more general formulae for 
derivatives of isotropic tensor functions derived by Chadwick and Ogden [5] . 

It is interesting to note that the tensor gradient function (first derivative) involves 
finite differences of the function / as well as its derivative evaluated at the eigenvalues of 
A. Similarly, the second derivative contains second order finite differences. The general 
result derived here shows that the coefficients for the n*'* derivative are related to an 
interpolating polynomial. 

Our main result is the following: 

Theorem 1.1. The n*'' derivative of the tensor function /(A) is given by 

d 

-V(")/(A)^ ^ f^:'l_^^^^M,^M,m...MAi^^,. (1.6) 

11,12, ...,i„ + i = l 

The coefficients /j'"]^, j are unaltered under all permutations of the n+l indices. The 

(n + 2)(n + 3)/2 distinct coefficients can be classified intc0 [-^^diiljhlj expressions f'^j^''^'' 
where Vi, Vj, i^k, are the number of occurrences of distinct indices i, j,k respectively, with 
i ^ j ^ k ^ i and Vi + Vj + = n + 1. The coefficient is 



f 



1 d-'-i fix) 

^ (Vl-\)\AX-^-^ n {x-araY^ 



l=i,j,k 



(1.7) 



Alternatively, the coefficient can be found from the unique interpolating polynomial P{x) 
of degree n that fits the data at the three points x = {a^, aj, ak} defined by the n + l 



^ [x\ is the floor function. 
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values f'^^Kai), f'-'^Haj), and f'^^^ak) for < / < i 
where /^'•'(a;) is the derivative. Let P{x) = PnX^' 

Jijk — Pn- 



-1,0 < J < i/j 



f 



The first few expressions for the coefficients are 
1 



ijk 



J ijk 



I ijk 



j.l,l.n-l 
I ijk 



= -J^''\o.k) 



1 



ak) 
1 



1=0 



E 



[n-l] 



(a, - akYf^'Hak) + («, - ak)f'{a,) 



1,0< K <i^k-l 
+ Po, tlicn 

(1.8) 

(1.9a) 
(1.9b) 

nf{aj] 
(1.9c) 



(ai - aj) \ {ai - akY 



ak) 



ri-2 



' ijk 



1=0 



{a, - Qffe)" 



-i-i 



{aj - akY'- 



-i-i 



(1.9d) 



(tti - Qffc)" 2 
(ttj - ttj) 



i^j 



ak)''' 
1 



1 + (n - 2) 



afe)"-2 

n— 3 ^ 



(ttj - ak) 



1 



-2-; 



[l + {n-l-2) 



(ai - Ckfc)" 
ai) 



'2-1 



a,- 



(1.9e) 



These are sufficient to determine all derivatives up to and including the fourth. Thus, 



the gradient tensor [n — 1) requires the two expressions /j^'^'^ and f^k^ evident from 



tO,lS 



(|1.4p : the second derivative (n = 2) requires Z^'^''^, fijl'^ and fljl'^, which may be read 
off from (|1.5P : the third derivative (n = 3) involves four distinct formulas for 



f0,2,2 , _fl,1.2 

fiik and - 



0,0,4 ^0,1,3 
ijk ^ J ijk ' 
0,0,5 



I ijk 



and the fourth derivative (n = 4) requires the five expressions /, 



ijk 



fijk^^ /yfe'"^' fv'k^ ^'^d Note that (HSg and ((Tjil) reduce to in the hmit 

as ai ^ ak- Similarly, (|1.9d|) and (|1.9ep reduce to (ll.9cp in the limit. 

The main objective of this paper is to prove Theorem ll.il Some new results concerning 
the properties of the fourth order gradient tensor V/ are also presented and the inverse 
fourth order tensor V~^/ is introduced. Both the gradient and its inverse are discussed 
with application to strain measure functions [HI [H] ■ The proof of Theorem 11.11 begins 
with a new derivation of the well known expression for the gradient V/. The essential 
structure of the second and higher order derivatives is shown to depend on a general 
algebraic identity. This identity also reveals the appearance of the characteristic finite 
difference terms. The proof is completed by making connections with contour integrals 
and with interpolation polynomials. The results are presented in terms of Kronecker 
products of tensors which makes the expressions more transparent. 
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The paper is laid out as follows. Notation is introduced in Section [2] followed by the 
derivation of Theorem II. II in Section [3l The inverse gradient tensor is introduced and its 
properties discussed in Section [H 

2. Notation and preliminaries. We consider second order tensors acting on vectors 
in a three dimensional inner product space, x — > Ax with transpose A* such that y- Ax = 
X • A*y- Spaces of symmetric and skew-symmetric tensors are distinguished, Lin — Sym 
® Skw where A e Sym (Skw) iff A* = A (A* = -A). Products AB GLin are defined 
by y • ABx = (A*y) • Bx. 

Psym is the space of positive definite second order tensors. Functions of a symmetric 
tensor can be phrased in terms of its spectral decomposition (jl.f p where A^ GPsym and 
the distinct eigenvalues a,, i = 1 . . . , c? < 3 are real numbers. The single function /(A) 

3 

is a special case of isotropic tensor functions of the form T(A) = ^ ji{a\^ai^a-i)Ki, 

i=l 

involving three functions of three variables. Chadwick and Ogden [2] derived first and 
second derivatives for more general situation (see [3] for a recent overview). Our interest 
here is strictly limited to one function / as in (|l.ip . 

Lin is the space of fourth order tensors acting on Lin. The square tensor product or 
Kronecker product ABB, LinxLin^ Lin, is defined in the usual manner as [llj 

(AHB)X = AXB*, VXe Lin. (2.1) 

The following property of H will be used extensively, 

(AKB)(XKY) = (AX)K(BY). (2.2) 

The generalization of (|2.1|) is 

( A^B . .M. . . C ) : XY^. . . Z = AXB*Y . . . ZC*. (2.3) 

n+l n 

Basic properties of the Kronecker product can be found in e.g. 

The gradient of a tensor function /(A) is a fourth order tensor V/ G Lin defined by 

V/(A) X = lim - [/(A + eX) - /(A)] . (2.4) 

Higher order derivatives are defined recursively in accordance with (|1.3[) . 

V("V(A)X= limi[v("-iV(A + eX)- V("-1V(A)], n>2. (2.5) 

We note some properties of the tensor gradient. First, the tensor gradient of the 
product of two functions is 

V[/(A).g(A)] = (lK,g(A))V/(A) + (/(A)HI) V.g(A). (2.6) 

For example, let g = 1/f = /^^, then using W[f{A)g{A)] — yields 

V/-i(A) = -[f-\A)^f-\A)]S/f{A). (2.7) 

The composition of two functions is / o g{A) ~ /(g(A)) has gradient 

V/o.9(A)=V/(g(A))Vg(A). (2.8) 
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Only symmetric tensor arguments are considered in this paper. The extension of the 
gradient tensor to non-symmetric tensor arguments is discussed by Itskov [7]. 



3. Proof of Theorem 11.11 The proof has several stages, starting with consideration 
of / in the form of a power series. This allows us to identify certain properties of the 
derivatives that are subsequently generalized to arbitrary /. Implicit in this approach is 
the idea that a function of a single variable can be approximated to any desired level of 
precision by a polynomial. The use of polynomials allows us to see the structure of the 
higher order derivatives, not only the first or second. 

3.1. A power series. For the moment assume 

f{x)^x"^, (3.1) 

then /(A + X) is 

(A + X)™ =(A + X)(A + X)(. . .)(A + X) 

=A™ + (A'"-iX + A"-2XA + . . . + XA™-i) 

-f- (A^-^X^ -i- A"'-3XAX + . . . + X^A"-^) + . . . . (3.2) 

The 0(X) and O(X^) contributions contain m and m{m — l)/2 separate terms, respec- 
tively. The term of 0(X''') has (™) elements. A more suggestive and useful form of (|3.2p 
is apparent by rewriting it as 



(A + X)" - A™ = ^ A™-'=XA'^-i + J2 A^^^^'^X A'-^XA'^-' + . . . 

fc=i fe=i 1=1 

m m—lk 

= Y^ (A™-'=KA'=-i)X+ ^ ^A^^i^'^HA'-^KA'^-' : XX 



k=l k=l 1=1 



m — n 



't i + ■ ■ ■ + * Ti + 1 — "m, — 

:XX...X + .... (3.3) 

n 

This permits us to consider the derivatives in sequence. Our objective is to understand 
the term of 0(X^^), but it is easier and more instructive to begin with the lowest, n = 1. 
3.2. The first derivative. Let us first focus on the 0(X) term. Consider the identity 

m m 

(AKI-IKA) ^A^^'^SA'^"^ = ^(A'"+^"'=HA''"^ - A™"''HA'=) 

k=l k=l 

= A"KI-IHA". (3.4) 

Then 

{Am - mA) [(A + X)'" - A™] = (A'"^! - IHA™) X + O(X^), (3.5) 
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which may be generaHzed to power series functions 

M 

fix) = (3-6) 

m=0 

as 

[Am - mA) [/(A + X) - /(A)] = [f{A)m ~ ih/(a)]x + o{x^). (3.7) 

Equation (|3.7p impUes 

Lemma 3.1. The tensor gradient of any differentiable function satisfies 

(AMI - IMA) V/(A) = f{A)mi - mf{A). (3.8) 

We will return to the form of this equation later when the inverse gradient tensor is 
discussed. At this stage we note an isomorphism between the tensorial identity (|3.4p 
and an algebraic equation in two variables. Let x and y represent the first and second 
occurrence of A in A""'' HA''"^ a;™-'=?/'=-i. Then eq. dO]) is equivalent to 



{x-y)J2 x"'''"i/'^ = ix- y){x"'-^ + + ^^'"-3^2 + , , , + j^m-i) 



fc=i 



a;'"-y". (3.9) 



This is the well known factorization of x™ — y™. The analogous relation for the second 
derivative of / is derived next. 

3.3. The second derivative. The algebraic equality analogous to dSll) for the 0(X2) 
terms involves three variables, say x, y and z, and is 

m— 1 k 

{x-y){y-z){z-x) x^-^- V"'^'"' = {z^y)x'^ + {x~ z)y'^ + {y~x)z"' . (3.10) 

k=l 1 = 1 

As proof of this statement note that the tensorial equivalent of (|3.10p is 

(AHISI - IKAHI)(IKAHI - IKIKA)(IHIHA - AHIKI)iv(2)/(A) 

= /(A)H(IKA - AMI) + (AH/(A)KI - IK/(A)HA) + (IMA - AKI)K/(A). 

(3.11) 

It may be checked, using p.3p . that this indeed holds. 

The algebraic identities for the n ~ I and n = 2 derivatives are, respectively, 

E-""V-^ = ^ + ^, (3.12a) 
^-^ X — y y — X 

^ m m m 

V Vx""i-V-^^'^'-' = 7 . + 7 \j r + 7 r. (3.12b) 

{x~y)[x~z) [y~x)(y-z) [z - x){z - y) 

This form of the identity will prove more useful in the general case of the n*'' derivative, 
considered next. 
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3.4. Higher derivatives. Referring to eq. p.3p it is clear that the term of 0(X") 
involves (™) products of the form 

A^'XA'^X ... XA*"XA'"+i, withi, + i2 + ... + in+i=m-n. (3.13) 

The analogous algebraic quantity has a closed form expression similar to that found for 
the first and second derivatives in (|3.12p . It is given by the following general identity 
among n + 1 independent variables: 

Lemma 3.2. 

m—n n + 1 „ 

53 X- . . . = . (3.14) 

H,i2,..^,i„ + l=0 1=1 Y{ {X, - Xj) 

i,+Z2 + ...+4„+i=m-ra 

As prooQ, consider the identity 



1 1 1 1 

{1 - Xlt){l - X2t){. ..){!- Xn+lt) ' ^ {1 - X^t) ' ^' ^ 



3 = 1 



The right hand side is simply the expansion of the left in partial fractions. Equation 
p.l4p follows by comparing the coefhcients of t'"^" on either side of (|3.15p . 

Motivated by the expansion p.3p . consider the following ansatz for the derivatives, 



V/(A) = /«A,HA„ (3.16a) 

d 

-V(2)/(A)= /^fcA,HA,KAfe, (3.16b) 

-V(")/(A)= Y /nl....,v+iA,,KA,,H...KA,„^,. (3.16c) 

...,i„,z„+i = l 

The key feature is the use of the spectral basis of A. 

The tensorial equations (|3.8p and (|3.1ip . or their algebraic versions p.9p and (|3.10p . 
imply 

(a, - aj)/5^ = /(a.) - /(aj), (3.17a) 

('21 

(ui - aj){aj - ak){ak - ai)/y-fe = ("fe - aj)/(aj) + ("i - "fe)/(aj) + ("j ^ ai)f{ak), 

(3.17b) 



■^Thanks to Doron Zeilberger for the proof. 
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These may be expressed in the obvious alternate form, 
(1) „ /(gp - fjaj) 



k 

J-(2) ^ 
J ijk 



i 7^ i, 



{ai - aj){ai - ak) {a^ - ak){aj - a^) (afe - ai){ak - aj) ' 



(3.18a) 
(3.18b) 

These relations do not, however, provide all the coefficients required to define the deriva- 
tives according to p.l6ap and (|3.16bp . Specifically, the coefficients f-p. 
are not defined, but may be obtained from p.l7p by taking the limits aj 
using I'Hopital's rule. The result is 



^ iik 



IP — f'{cx.i) in agreement with \1A\ 

fja,) /(a.)-/(a,) ■ , 

Cti-Ctk (Cti— Ofc)^ ' ' ' 

i/"(a,), I = k, 



f^^^ and f^^^ 
hkk ^^'^ Ikkk 

-f Ui etc. and 
and 



(3.19) 



in agreement with Theorem 11.11 and with (jl.Sp from [Tir . 

How does this generalize to the n*'* derivative, and is it correct to take the limits as 
described to find the "missing" coefficients? In partial answer to the first question, note 
that Lemma [3.21 applied to the fmiction / represented as a power series, p. 61) . suggests 
the equivalence 



E 



n+l 



An) 



n > 3, 



(3.20) 



where the Xi and Xj are replaced by ai and aj. However, no matter what combination of 
indices are taken for the coefhcients on the right hand side, there will be pairs of indices 

(3) 

with the same values, e.g. fijl^- In this sense the cases n = 1 and n = 2 are special 
because it is possible to have indices that are distinct, which enables us to write the 
identities (|3.18p . This does not apply for n > 3. The indicial overlap must be taken into 
account in order to arrive at the correct form of (|3.20p . Not surprisingly, the correct 
result is equivalent to taking the appropriate limit of the left member of p.20p . Consider 

(3) 

the third order derivative, for instance. The coefficient corresponds to xi = at, 



X2 



Qffe in the limit as x^ 



ak- That is. 



ak){a.j - x){aj 



ai) 



{ui - aj){ai ~ ak){ai ~ x) (< 

fio^ ^ 

(afc - ai){ak - aj){ak - x) {x - ai){x - aj){x - ak) 

' f{ai) - ,f{ak) /(aj)-/(afe)] ^ /'(a/c) 



(a, - aj) 



(a^ - Qffc)^ 



Oik) 



{ai - ak){aj - ak) 



(3.21) 



in agreement with (jl.9dp for n = 3. The remaining coefficients fil\k^ fukk ^^'^ fkkkk 
be obtained by taking further limits. This method of evaluating the coefficients, while 
correct, is particularly tedious as n becomes larger and more limits must be taken. 



f(3) 



can 



Let us therefore consider the form of /, 



(n) 



' «1 ,l2,---,'ln + l 

is unchanged under any permutation of the indices ii, 



. First, it is clear that the coefficient 



«2, ■ 



In other words, it is 
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totally symmetric in the n + 1 indices. Secondly, the indices assume only three different 
values: 1, 2 and 3. Since there are n + 1 indices then apart from the first and second 
derivatives fj;p (n = 1) and f^^^ {n = 2) there will always be coincident indices present 

in] 

in the coefficient f. , , . These are indices with the same value whether that value 
is 1, 2 or 3. For n > 3 there are at least n — 1 coincident indices. Suppose the index k, 
distinct from i and j, occurs > 1 times in {ii,i2, . . . ,in+i\- Similarly, the indices i 
and j each occur Vi> Q and Vj > times. This suggests the change of notation 

itl,....„+, - with + + n + (3.22) 

We are now ready to describe two approaches for calculating these coefficients. 

3.5. A closed form expression. We first note that the coefficient /i"]^ ^ of (I3.16cp 
can be evaluated as a contour integral. This follows from (|3.20p and the integral identity 

E"^^ f(xi) 1 f dzf(z) /„ „„x 

„+i = ^ y „+i (3-23) 

n C n (Z-Xl) 

j=l 1=1 

where C is any contour in the complex plane that encloses the set of points {xi, X2, ■ ■ ■ , Xn+i}- 
The integrand in (j3.23p must be modified to account for multiple occurrences of the in- 
dices, in the notation of (|3.22p . We obtain the integral identity 

r'^^,l^,,l^k 1 f dz/(z) 



' 27ri J {z- aiY^ (z - UjY^ (z - auY" ' ^^'^'^^ 
c 

which can be evaluated by residues, with the final result l|1.7p . 
3.6. Interpolation interpretation. The following result is key: 



Lemma 3.3. If {xi, X2, . . . , Xn+i\ are n+1 distinct points then 

n\ dx 



y Hxj) ^ i.^P(:,), (3.25) 



i=l 



n {Xi - Xk) 



k = l 

where P{x) is the Lagrange interpolating polynomial of the n+1 function values f{xi), 
i = 1,2, . . . ,n + 1. 

The proof follows from the explicit form of the interpolating polynomial of degree n, 

n+1 

2=1 

where Li(a;), i = 1, 2, . . . , n + 1, are the Lagrange polynomials satisfying Li{xj) — Sij, 

n+1 I n+1 

U{x)^Yi^x-xi)r^{x,~Xk). (3.27) 

3 = 1 ' fc=l 

Differentiating P{x) of ([H?^ yields the identity ([323) • 
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This connection between the coefficients and polynomial interpolation provides an 
alternative means to calculate the coefficient /j^tJ ■ The interpolation problem that 
must be solved is not the straightforward one of Lemma 13.31 with a single datum at 
each of n + 1 distinct points. Some of the points can have more than one piece of data. 
Specifically, associated with x = ai where I = i, j and fc, will be the i>i data comprising 
f{ai), f'{ai), . . . , /('''"I) (a;). In this way there is a well defined and unique interpolating 
polynomial of degree for each f^J^!^' with 

The interpolation problem that must be solved does not fit into the standard category 
of Lagrange interpolating polynomials, since some or many of the points are coincident. 
However, the procedure for generating the polynomial is not difficult. Appendix |A] il- 
lustrates the method by example, for the case of / known at three points along with all 
derivatives at one of these points up to 

3.7. The number of independent coefficients. The final detail is a count of the number 
N{n) of independent expressions of the form f^j!^^'^'° required for the v}^ derivative. 
This is the number of elements in the set of integer triples {/, J,Ke Z: 0</<J< 
K, I + J + K = n + l}. That is, 

N{n) = dim {/, J e Z: < I < J < n + 1 - I - J}. (3.29) 

This defines a region A{n) on the integer grid of {/, J} bounded by the three straight 
lines Li : I — 0, L2 I — J = and L3 : I + 2J = n + 1. N{n) is the number of points 
in A(n), including on its boundary. Enumeration yields 

»(„) = l(!l±if±ij (3.30, 

where the floor function [a;J is the integral part of x. 

This completes the proof of Theorem 11.11 Note that the two methods discussed for 
calculating the coefficients, based on the contour integral and the interpolation problem, 
are both related to divided differences [4]. 



4. The inverse gradient. 

4.1. Definition. The inverse tensor function V^^/( A) G Lin is defined by (V/)V^^/ = 
(V^^/)V/ = I. As an example application let /'"^ be the inverse of /, such that 
/ o /''^^(A) = /'"^ o /(A) = A. It follows from dSiH]) that 

Vr^(A)-V-V(r^(A)). (4.1) 

It is clear from (|1.4I1 that the definition of V^^/(A) is problematic if f'{ai) vanishes. 
This possibility is precluded by restricting consideration to strictly monotonic strain 
measure functions acting on positive definite tensors. 



HIGHER DERIVATIVES AND THE INVERSE DERIVATIVE OF A TENSOR- VALUED FUNCTION OF A TENSOR 



4.2. Application to strain measure functions. The function / is a strain measure [6lll2j 
if it is a smooth function / : K"*" M which satisfies 

/(1) = 0, /'(1) = 1, /'>0. (4.2) 

For the remainder of the paper we restrict attention to strain measure functions. Exam- 
ples win be presented for the Seth-Hill strain measure functions, 

f^"'\x) ^m-\x"' -I). (4.3) 

The inverse /'"^ is well defined since strain measure functions are one-to-one. Hence, the 
inverse gradient is relevant to calculating the gradient of the inverse function, via the 
identity f^ . 

Strain measure functions act on positive definite strain tensors S G Psym, with 

d d 

S = ^A,S„ 1 = 51^*' •^*>0' S,S, -SjS, = <5,jS„ (4.4) 

1=1 1=1 

where d < 3 is now the eigen- index of S. 

Accordingly, we extend the concept of positive definiteness to fourth order tensors. 
It is necessary to define an inner product on Lin, which is done in the usual manner 
as A • B — tr(AB*). Then Lin is the space of fourth order tensors acting on Lin, 
X AX with transpose A* such that Y • AX = X • A*Y for all X, Y eLin. The vector 
space may be decomposed Lin = Sym© Skw where Sym and Skw denote the spaces 
of symmetric ( A* = A) and skew-symmetric (A* = — A) tensors, respectively. Any 
A G Lin can be uniquely partitioned into symmetric and skew parts: A= aW + a(-), 
where A^=^^ = ( A± A*)/2. The identity I satisfies IX = X for aU X e Lin. The product 
AB e Lin is defined by Y • ABX = ( A*Y) • BX. Psym is the space of positive definite 
fourth order tensors: A S Psym iff X • AX > 0, for all nonzero X G Sym. 

The spectral form of A e Psym is 



D D 

^aiAi, = 
1=1 1=1 



a/ A/, I=>^A/, A/Aj=<J"' ^ (4.5) 

I J, 



where D is the eigen-index, a/ > and A/ G Psym. Functions of fourth order tensors 
can therefore be defined according to 



D 



/(A) = > /(a/)A,. (4.6) 



1=1 



For any A ePsym and its eigen-tensors A^, i = 1, . . . , d of (|1.2p . define the associated 
set oi D = ^d{d + 1) basis tensors for Sym by 

(4.7) 

,D = ld{d+l). 

Here / = d + 1, . . . , D corresponds to distinct pairs (i, j) with i < j. Thus, _D = 6 for 
d = 3, while D — 3 for d — 2 and D = 1 in the trivial case of c? = 1. It may be readily 
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checked that A/ Aj = AjSij. The identity I = I Kl I imphes the partition of unity 

1=1 

In particular, 

Lemma 4.1. The gradient of a strain measure function and its inverse are positive definite 
fourth order tensors, i.e. V/(A), V^^/(A) e Psym, and 

(V/)V-V = (V-V)V/ = I. (4.8) 

The spectral forms are 

D D 

V/(A)=^/,A,, V'V(A) = ^/7^A,, (4.9a) 
1=1 1=1 

with 



/'(«/), /=l,...,rf, 

/ = rf+l,...,D. 



The proof is evident from the spectral decompositions. The positivity of //, / = 
1, . . . , D is a consequence of > 0, i — 1, . . . , d combined with the monotonicity of / 
and the mean value theorem. The positive definite nature of V/(A) implies that it has 
a unique positive definite inverse. 

The gradient function and its inverse has an interesting alternative representation for 
integer and fractional values of m, respectively: 

Lemma 4.2. For integer values of m 7^ 0, 



J2 A"-'=KA'^-\ m>0, 
V/["l(A) = ^in (4.10a) 
\m\ j2 A"'-''^A^-\ m<0. 

[^k=m+l 

EA1--HA---, m>0, 

. . . (4-lOb) 
E A^-^SA^-^, m<0. 

The proof uses the expansion of eq. (|4.9p with // ^ /j™' where 



r„,i ( aV* ^ , / = 1, . . . , d 
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Consider the case of positive m. Then using the spectral form for A 



... ~k k~l 

a,- a,- 

m ^ — ' m ^ — ' " ^ — ' 

k=l = l k=l 



d 

a, N k 



2—1 ij = l 



TO ^ — ' ■' '■ — ' a," 

:.3 = i k=\ 



= V/['"1(A), (4.12) 

where the identity a; + + . . . + x™ = x(l — a;™)/(l — x) has been used. The proof for 
negative to is similar. The results for V^^/'^™^(A) are a consequence of the identity 

/j-\A) = l(A^), TO^O. (4.13) 
Some examples of interest, first for integer values of to: 

V/(i)(A) = IKI = I, (4.14a) 
V/("i)(A) = A"iHA"\ (4.14b) 

V/(2)(A) = -(IHA + AKI), (4.14c) 



V/("^)(A) = i(A"2g]A-i _^ A"^SA"2-)_ (4.14d) 



The inverse tensors 



V^V'^HA) = I, (4.15a) 

V^V^^^^A) = AHA, (4.15b) 

follow by observation, as do the general relations 

V/(-™)(A) = A-"HA-" V/["1(A), (4.16a) 

V-ij(-m)(A) = A'^HA" V-V[™1(A), (4.16b) 

which express tensors in terms of their counterparts of opposite sign. 

AX + XA = C is an important and common equation in mechanics |131 lllj . The 

solution may be written simply and succinctly as X = (AKII + IK1A)^^C where [S] (eq. 
(22)) 

d 

(AHI + IHA)-i y^ {a,+oij)-^K.mAj. (4.17) 

ij = l 

In the current notation this may be expressed X — ^V^^/'^' (A)C, see eq. (|4.14cp . In 
fact, the identity ()4.10ap implies that the solution to 

A'"-iX + A^-^XA + A'^-^XA^ + . . . + XA™-i = C, (4.18) 
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IS 



X ^ lv-V['"l(A) C = ^ (^^^)A,^A, C. (4.19) 
m ^-^ a" — a" 

The solution to AX+XA = C can be written in a form that does not require the spectral 
split, e.g. X = [2(/i/2 - h)]-^ [{If - I2)C + ACA - (A^C + CA^) + hl^A-^CA-^] , 
[TT] (eq. (27)) where Ii,l2,h are the invariants of A. Similar expressions, although far 
more complicated, can be generated for the solution to (|4.18|) and related equations, but 
we leave such matters aside. 

Examples of fractional powers include the well known case of i llj, for which (j4.10bp 
gives 



\7-^fih)[A)^]^{^^Am + m^^X) => V\/A = (VaKI + IH\/A)-\ (4.20) 
imilarly, 

V-V(-3)(A) = i(\/AKA + AS\/A) VA"i/2 ^ -(\/AKA + AE\/A)-\ (4.21) 



and 

y-i/^^^A) = i(A^/^HI + A^/^KA^/^ +IHA^/^), (4.22a) 

y-i/^^^^A) = i(Ai/3HA + A^/^HA^/^ + AHA^/^), (4.22b) 

imply expressions for VA^^^"^, etc. 

Note that the case of m = 0, corresponding to 

V/(*')(A) = Vln(A), (4.23) 

is specifically excluded in Lemma 13.31 Although there is no expression for V ln(A) 
analogous to (I4.10ap . its inverse can be represented as follows: 

Lemma 4.3. The tensor V^^ ln(A) — V^^/("^(A), which is the inverse of Vln(A), is 

1 

V"^ ln(A) ^ j dxA'^MA}-'', for any A e Psym. (4.24) 



The proof follows from (j4.10bp by taking the limit of to oo. An alternative and 
direct proof can be obtained using the spectral form of A, 

1 d 1 

j dccA^HA^-^ = J2 A.,KAj- J dxa^a^ 



Q!j — a 
In Ui — In q; 



E I7^EfcA.^A„ (4.25) 



where the ratio becomes if i = j. The right member of (j4.25p is obviously the inverse 
of Vln(A) from eq. (fTi]) . 
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4.3. The tensors J, J*, IK, and the solution to AX — XA = Y . Consider the equation 

AX - XA = Y, A e Psym, (4.26) 

for the unknown X in terms of Y which is either symmetric or skew, and X is of the 
opposite parity [5] . The equation can be written 

J(A)X = Y, (4.27) 

where 

JI(A) = AHI- IHA, A e Psym. (4.28) 

We will only consider J(A) for symmetric A, implying J G §ym and J maps Sym^Skw 
and Skw— +Sym. Therefore, J does not possess eigenvalues, eigenvectors or an inverse in 
the usual sense. However, it is possible to define the pseudo-inverse, or equivalently the 
Moore-Penrose inverse, J* , which satisfies 

JJJJ* J = J, 

rir = r. (4.29) 

Further understanding comes from consideration of the spectral form 

d 

J(A) = ^ (a. - aj)AMAj- (4.30) 

ij = l 

This representation implies that A^ are null vectors, i.e. J(A)Ai — 0, i — 1, . . . , d. 
Noting that 

J(A)[AjKlAj ± AjKlA,] = (a,; - a.j)[AMJ^j T AjKlAj], (4.31) 

gives 

r(A)[A,KlAj ± AjKlA,] = {a^ - aj)-^[AMAj T AjKlA,], for i ^ j only. (4.32) 

The caveat i ^ j crucial, and indicates that the non-null parts of J and its pseudo- 
inverse define maps between three dimensional subspaces of Sym and Skw. In particular, 

d 

r (A) = (AHI - IB A)* = ^ (a, - aj)-^AMAj. (4.33) 

This formula for the pseudo-inverse clearly satisfies (|4.29p . It is interesting to compare 
the form of with (jiTZ)) . 

Note that 

d D 

J*J= jr = I-^ A/ = Aj. (4.34) 

1=1 I=d+1 

Define K(A), A ePsym, and its pseudoinverse, 

d d 

]K(A)=^a,A„ ]K*(A) = ^ar^ A, = K(A"i), (4.35) 

i=l i=l 

then 

J(A) r (A) -I- K(A) K* (A) = I for any A e Psym. (4.36) 
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The tensors Jf and K allow us to write V/(A) and its inverse in a succinct manner as, 
respectively, 

V/(A) = IK(/'(A)) + r(A)j(/(A)), (4.37a) 
V- V(A) = K* (/'(A)) + J(A) r (/(A)) . (4.37b) 



Appendix. 



A. Solving the interpolation problems. Consider three points, say xi, X2 and X3. 
Given f{xi), /(X2) and /^'^(xa), I = 0,1,2, ... ,n — 2 find the interpolating polynomial 
P{x) of degree n. For simplicity, but no lack of generality, take 0:3 = 0, and consider the 
ansatz 



P{x) 



1=0 



Pi 



1^1,2,. 



2. 



(A.l) 



This satisfies all n — 1 conditions at a; = 0, and the conditions at xi and X2 are met if 
Pixi) = f{xi) and P{x2) = f{x2)- These imply a pair of linear equations in and 
Pn which are easily solved. Only the coefficient p„ is required as this determines the n*^ 
derivative of P{x), 



1 



Pn 



Xi - X2 



n-2 

E 



n {x" 



(A.2) 



This implies eq. (jl.9dp . 

Consider the case where the values of the function and its first n — 3 derivatives are 
given at one point, say x = 0, along with f{xi), f{x2) and f'{x2). Using the same 
ansatz as (|A.ip with the three unknowns Pn-2, Pn-i and p„, implies three simultaneous 
equations. The solution for p„ is 



1 



Pn 



{Xl - X2)^ 

n— 3 



E 

1=0 



n-2 

n 



X2 



„n-2 



„n-2 



—^{X2~xi) 



[l-(^-2-0(--l)] 



X2 



(A.3) 



This yields the expression (|1.9ep . 



Acknowledgment. Thanks to Doron Zeilberger for being true to his motto: "Who 
you gonna call?" . 
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